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Abstract 

We have derived the quark mass dependence of m 2 , ttiawi an d f-^, using the chiral perturbation 
theory which includes the a 2 effect associated with the explicit chiral symmetry breaking of the 
Wilson-type fermions, in the case of the Nf = 2 degenerate quarks. Distinct features of the 
results are (1) the additive renormalization for the mass parameter m q in the Lagrangian, (2) 0(a) 
corrections to the chiral log (m q logm q ) term, (3) the existence of more singular term, log m q , 
generated by a 2 contributions, and (4) the existence of both m q \ogm q and logm g terms in the 
quark mass from the axial Ward-Takahashi identity, th-awi • By fitting the mass dependence of m 2 
and mAwi, obtained by the CP-PACS collaboration for Nf = 2 full QCD simulations, we have 
found that the data are consistently described by the derived formulae. Resumming the most 
singular terms logm g , we have also derived the modified formulae, which show a better control 
over the next-to-leading order correction. 
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I. INTRODUCTION 



One of the most serious systematic uncertainties in the current lattice QCD simulations is 
caused by the chiral extrapolation. Due to the limitation of the current computational power, 
one can not perform simulations directly at the physical light quark(up and down) mass. 
Instead, one has performed simulations at several heavier quark masses and has extrapolated 
results to the physical quark mass point, using the polynomialflinear, quadratic, etc.) or the 



formula derived from the chiral perturbation theory (ChPT)jJJ. These extrapolations cause 
large systematic uncertainties, in particular in the case of full QCD simulations, where the 
lightest quark mass employed in the current QCD simulations is roughly half of the physical 
strange quark mass (m n /m p ~ 0.6). 

Recently more serious problem has been pointed out, in particular for full QCD simula- 
tions with Wilson-type quarks: the expected chiral behaviour predicted by the ChPT has 
not been observed. For example, the behaviour of the pion mass m 2 as a function of quark 
mass m q is given by 

ml = Am q [l + l6 ^ /2 log(Am q /A% (1) 

where A is some scale parameter. Since the pion decay constant is experimentally known as 
f n = 93 MeV, only A and A are unknown parameters. Unfortunately, such a two parameter 
fit can not explain lattice data well, which looks almost linear in the simulated range of 
quark masses. If one includes /„- as a free parameter, the best fit typically gives f% > 5x 
(93 MeV) 2 Q. 

The most widely accepted interpretation for this discrepancy is that the simulated range 
of quark masses in the current simulations is still too heavy to apply the ChPT. If this 
interpretation is true, the current lattice simulations with the (Wilson-type) dynamical 
quarks lose a large part of their powers to predict properties of hadrons at the physical light 
quark masses. 

In this paper, we investigate a theoretically more natural alternative that the explicit 
breaking of the chiral symmetry by the Wilson-type quark action modifies the formulae 
of the ChPT at the finite lattice spacing. We first derive formulae in the modified chiral 
perturbation theory for the Wilson-type quark action, denoted by WChPT in this paper. 
Such attempts have been made before at the leading order jjj and the next-to-leading order 0. 
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At the leading order |3|, the WChPT predicts the existence of the parity-flavor breaking phase 

n h q 

transition^, lalU f° r the 2 flavor QCD as long as massless pions appear at the critical quark 
mass. This analysis has also shown that the 0(a 2 ) chiral breaking term play an essential 
role to generate the parity-flavor breaking phase transition, which is necessary to explain 
the existence of the massless pions . In the next- 

to-leading order analysis^, however, only the 0(a) breaking effects are included, and it is 
concluded that the effect of the chiral symmetry breaking can always be absorbed in the 
redefinition of the quark mass, so that all formulae in the ChPT remain the same if one 
replaces the quark mass m q with m q — m c , where m c is the additive 0(a) counter-term for 
the quark mass. In the section |H1 we perform the next-to-leading order calculation in the 
WChPT including 0(a 2 ) chiral symmetry breaking effects. To make the difference between 
WChPT and ChPT clear, we consider only the case of the Nf = 2 QCD with degenerate 
quark masses, and derive the formulae for mass and decay constant of the pion as well as the 
axial Ward-Takahashi identity quark mass, as a function of the "quark mass" in the effective 
theory. In section the derived formulae are applied to data of pion mass and the axial 
Ward-Takahashi identity quark mass calculated by the CP-PACS collaboration 0. We show 
that data are consistent with the formulae. We have attempted the resummation of the 
most singular term, and have derived the modified formulae in section IIVI Our conclusions 
and discussions are given in section IVl 

II. WILSON CHIRAL PERTURBATION THEORY 
A. Derivation of effective Lagrangian 

It is difficult to derive the effective chiral Lagrangian for mesons directly from lattice 
QCD with the Wilson-type quarks using the symmetry, since the quark mass requires a 
counter term m c , which diverges as g 2 /a near the continuum limit, so that m c a = 0(1) 
and the conventional power counting of a fails. Therefore, following the proposal^, [^J, we 
overcome this problem by first matching the lattice QCD to an effective continuum-like QCD 
including the scaling violations into higher dimensional local operators ^, then match the 
latter to the effective Lagrangian for the Wilson chiral perturbation theory(WChPT). 

Close to the continuum limit, the lattice QCD can be described by an effective action in 
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the continuum, which is expanded in power of a as 

^cff = 5*0 + aSi + a 2 S 2 H , (2) 

where Si contains chiral non-invariant terms only, while S2 contains chiral invariant as well 
as chiral non-invariant terms. By using the equation of motion and the redefinition of the 
quark field, quark mass and the coupling constant, only one term is relevant in Si: 

Si = antpa^F^tp H . (3) 

The similar analysis can be done for So 

We now derive the effective Lagrangian of the WChPT from S e fr, using the symmetry 
of S e ff such as parity, axis inter-change symmetry (rotational invariance in the continuum 
limit), and the chiral symmetry. The last one is explicitly broken not only by the quark 
mass m but also by the breaking terms in Si and S2, whose coefficients are denoted as r»( 
i = 1, 2, 3, ■ ■ •). One can make S e ff formally chiral invariant by transforming m and r^'s to 
compensate the chiral variation of ip and ip. For example, if one writes the quark mass term 
as 

c.MP, ; c + cM'Plc. (4) 

this term is invariant under 

V -> (RPr + LP L ) z/>, ^ - i> (tfP R + rfP L ) (5) 
M -> l.MIV. M f P\Pi:. (6) 

where R and L are SU(iVj) chiral rotations. The usual mass term is recovered by setting 
M = = m. The similar transformations can be defined for r^'s, but we do not give them 
explicitly since the detail of them is irrelevant for later discussion. From this argument one 
concludes that the effective Lagrangian of the WChPT should have this (generalized) chiral 
SU (N f ) R ® SU(N f ) L symmetry. 

As mention in the introduction, we consider the Nf = 2 case to make our argument 
simple and clear. In this case, the chiral field for the pseudo-scalar mesons(pions) is given 
by 



E(a;) = £ exp |^7r a (x)t7/j = E [cos(tt//) + i7r a t a sin(7r//)] 



(7) 



where n a (x) is the pion field, t a = a a is the ordinary Pauli matrix and / is the pion decay 
constant, whose experimental value is 93 MeV. The norm and the unit vector of the pion 
fields are given by n 2 = n ■ n = J2a 7! ~ a ' n ' a , an d vr a = ^ a /^, respectively. As discussed in 
ref. jsj], the vacuum expectation value E may have a complicated structure, leading to the 
spontaneous breaking of parity-flavor symmetry, but in this paper, we stay in the phase 
without this symmetry breaking, so that Eo = 12x2- Under the chiral rotation, this field is 
transformed as E —>■ LUR' '. Under the transformation that 7r — > — n, called "parity" in this 
paper, 2^ — > Et. 

Using this field, we define the following naive operators for Scalar (S), Pseudo-scalar (P), 



Vector (V), and Axial- vector (A): 

S° = -tr (E + E^ = cos(tt//), S a = -tr t a (E + E f ) = (8) 

P° = -tr (E - E f ) =0, P a = ^tr t a (e - E f ) = ii a sin(7r//) (9) 

Lj = itr (E^Et) =0, L« = itr t a (E^Et) (10) 

i$ = ^tr (Et^E) =0, R; = itr t a (st^s) (11) 

K = \(i<l + K) = Q > A l = \( L l-rf) = Q ( 12 ) 

V; = \(Ll + R a ,) = 2e afec 7r & sin(vr//)^[7r c sin(vr//)] (13) 



A; = ^(L^-R;)= i {fr a sin(rc/f)d,[cos(n/f)]-cos(7r/f)d^ (14) 

where the suffices and a mean the flavor singlet and triplet, respectively We also introduce 
Left-handed(L) and Right-handed(R) currents for later use. Due to the speciality of the 
Nf = 2 case, some of the above operators are identically zero. Here we do not consider the 
Tensor (T) operator, which must contain two derivatives, since it does not contribute to the 
1-loop calculation in this paper. 

Now we construct the effective Lagrangian, which must be invariant under parity, axis- 
interchange symmetry and the (generalized) chiral symmetry. In the 1-loop calculation, 
which gives the main contribution at the next-to-leading order in the chiral perturbation 
theory, it is enough for us to construct the effective Lagrangian up to the order m, where 
m is the quark mass in the effective theory. On the other hand, we must include the 0(a 2 ) 
effect to realize the massless pions at a ^ [3]. At the next-to-leading order, 0(m 2 ) counter 
terms (Gasser-Leutwyler coefficients) are also needed. We do not include, however, these 
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terms in our effective Lagrangian, since we will not intend to determine them in this paper. 
Instead we introduce arbitrary scale parameters in log(m) terms which appear in 1-loop 
integrals. Roughly speaking, we consider the situation that l^>a>m~p 2 >a 2 > ma ~ 
p 2 a > m 2 ~ p 4 ~ mp 2 , so that all terms up to ma or p 2 a in this inequality will be included 
in the effective Lagrangian. 

The chirally invariant contribution at the leading order, which has the least number of 
derivatives, is constructed from L a or R a as follows: 



2 £ = 2 £ R;r; = tr [^E^E 

a=l a=l 

= 2{^[cos(7r//)]^[cos(7r//)] + ^[7r a sin(7r//)]^[vr a sin(7r//)]}, (15) 



= « = 0. (16) 

Note that R a ^L a ^ term is prohibited by the parity invariance. The chirally non-invariant 
parity-even term accompanied with one power of m, r\ = 0(a) or rj> 2 = 0(a 2 ) is uniquely 
given by S°. The chirally non- invariant terms whose coefficients include r\ = 0(a 2 ) or 
T\-m— 0(ma) are given by (S* ) 2 , J2a{P a ) 2 or tr(E + S''") 2 . For the Nf = 2 case, however, 
the latter two terms are not independent, as evident from the expressions that Da=o(-P a ) 2 = 
(/S ) 2 — 1 and tr(E + S''") 2 oc (S ) 2 . An independent term at 0(ap 2 ) is given uniquely by 
S x tr^E^E], since tr[(E + E^E^E] is not independent for SU(2). 

Gathering all terms up to m,p 2 , a 2 and ma,p 2 a, the effective Lagrangian becomes 

L cS = £ [l + c (S° - 1)] tr {d M E%E} - Cl 5° + c 2 (S ) 2 , (17) 

where parameters Co, ci and c 2 have the leading m and a dependences as 

Co = W a + O(m) (18) 
d = Wia + Bim (19) 
c 2 = W 2 a 2 + y 2 ma + 0(m 2 ). (20) 

Since Co is dimensionless and Cx and c 2 have the mass dimension 4, Wo ~ A(l + O(Aa)), 
Wi ~ A 5 (l + 0(Aa)), W 2 ~ A 6 (l + 0(Aa)), F 2 ~ A 4 (l + O(Aa)), B 1 ~ A 3 , where A 
represents some mass scale of the theory such as Aqcd- The (sub-leading) a dependence 
of these parameters comes from the chiral breaking terms of a 2 5* 2 in the effective action 
eq.(J2J), which correspond to rj> 2 = 0(a 2 ) terms in c and Cj, or r 1 ■ rj> 2 = 0(a 3 ) and 



in ■ Tj> 2 = 0(ma 2 ) terms in c 2 . Chirally invariant parameters such as / receive 0(a 2 ) 
corrections from chirally invariant 0(a 2 ) terms in a 2 S 2 - Note that W , Wi, V 2 ~ 0(a) if 
non-perturbatively 0(a,ma) improved fermions are employed for the lattice QCD action. 

For later use, we define the operators in the effective theory, which correspond to the 
ones in QCD up to 0(a): 



= Z S S°{1 + c s (S° - 1 )}, [P a ] = Z P P a {\ + c P (S° - 1)} 



V 



A a 



cvd„S°, 



V!; 



Z v V«{l + c v (S°-l)}, 



(21) 
(22) 

= Z A {A; (f + c A (S° - 1 )) + c A d,P a ) (23) 

where cs,p,v,a and c A ,v are 0(a) in general, or 0(a 2 ) if the lattice action and operators are 
non-perturbatively 0(a) and 0(ma) improved. 



B. Next-to- leading order calculations 



To perform the next-to-leading order (1-loop) calculation, we expand L e fr in terms of the 
pion field ir a as 



L cS = const. + - 



ci - 2c : 



P 



2 ^r 2 
7T 



+ 



6P 



(7r • d^n) 2 - (1 + -c )(df,TT ■ d^u)^ 



+ 



7T 



2\2 



4!/< 



lc 2 - Ci) 



and the operators as 



5° 



V 



4 a 



7T 



7T 



1 - 



7T 



2!/ 2 



■ 7 7T" 



1 - 



71 



3\f 2 



;i + 3c P ) 



abc " M 



7T 



/ 2 



3!/ 2 



;i + 3 Cy ) 



/ V 3!/V 
(1 + ^)- ^(l + — r -) + 3 (1 - - ) 



Using the pion propagator at the tree-level, which is given by 

1 



K(-PK(P)>0 = Sab 



P 2 + WIq 



m n = 



ci - 2c 2 



(24) 

(25) 
(26) 
(27) 
(28) 
(29) 

(30) 
(31) 
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we evaluate loop integrals as usual 
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Tti in 

(d^ a {x)d v 7t b {x)) = 5 a M (-mil) , 



(32) 
(33) 



where we introduce an arbitrary scale parameter A resulting after removals of power diver- 
gences of loop integrals by the local counter terms. Therefore, although we use the same 
symbol, this A varies depending on physical observables. 

The inverse pion propagator at the 1-loop level is calculated as 



r(2) _ 
^cS — 



6/ 



where 



(9„7r fl ) +m 7r 7r R 



7T = Z^ R 



z = 



1-^(2 + ^) 



ml = m 2 



3/ 



1 + 



mt 



;i + 6c )log-^ + 



32tt 2 / 2 

For the axial- vector current, we obtain 

iZ 
1 

therefore the decay constant at the 1-loop order becomes 



ml 5ca_ ml 
A 2 16tt 2 / 4 8 A 2 



m\(x)n b R (y)) = 5^(d,n R (x)n b R (y)) Z 1 / 2 [(1 +c a )-^(a + 



9c A - 3c A 



iZ , 



Taking Z A = —i^/2f 2 , we have 



1 - 



167T 2 / 2 



1 + 



3c A Hca 3cq 
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4 7 l0g A^ 



A = /(1 + ca) 



m„ . . , ml 

16# (1 + C '- )lDg A? 



(34) 

(35) 
(36) 

(37) 
,(38) 
(39) 

(40) 



where cy^ = 3c^/2 — llcU/6 — 3c /4. Note that receives an 0(a) correction even in the 
chiral limit: f n = /(l + £4). 
Similarly, we have 



(dM»>R(y)) = (ir R (x)ir b R (y)) V2fmlZ^ 

1 - 



d + ~c A )-^(4 + 9CA - 3 ~ CA 



([P a }(x)n b R (y)) = iy(v a R {x)K b R {y))Z^ 



5/ 
3!? 



3/ 2 
(1 + 3c P ) 



(41) 
(42) 
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Then the PCAC quark mass mAwi is given by 



<[p°](*K(y)> 



iZp 



ml(l + c A ) 



rrin 



32vr 2 / 2 



^1 + 30^ - 11ca/3 - 5c P ) log 



A 2 



l + c A 



-m r 



1 + 



m O C m A wi + 10c 2/ P 



log 



^0 

A 2 



32tt 2 / 2 

where 1/(2B ) = \f2~P/(iZ P ) and c mAWI = 6c - 3ca + llc^/3 + 5c P . 
Let us recall the leading m and a dependences of the parameters: 

c = W a, C\ = W\a + Bim, c 2 = W 2 a 2 + V 2 ma 
cp = Wpa, ca = Wacl, ca = Wao>, 

and then the pion mass at tree level is written as 

ci - 2c 2 m{B 1 - 2V 2 a) + aW x - 2a 2 W 2 



(43) 



(44) 
(45) 



m r 



where 



.4 



P 



B 1 - 2aV 2 



P 



A(m — m c ) = Am 



R 



(46) 



m. 



Wx - 2aW 2 



m R = m — m c . 



(47) 



P ' "" c B l -2aV 2 1 

Here it is noted that m c = 0(a) does not correspond to 1/(2K C ) in lattice QCD, since the 
1/a contribution to the quark mass is already subtracted in m. Furthermore, for m < m c , 



pion would become tachyonic (m 2 , < 0). As discussed in ref. j^, however, as long as Co = 

W 2 a 2 + V 2 m c a = 0(a 2 ) > o, the uu™ s „ wi g P w —nEaa 

occurs at m — m c = 0(a), so that m 2 , is always positive. In other words, the 0(a 2 ) 
contribution in c 2 is necessary for the consistency between the PCAC relation (m 2 ~ m q ) 



and the absence of tachyons 



We summarize the result of the 1-loop calculation in terms of m R and a: 

,2 



rriAwi 



P 



Am R 
A m R 

f(l + ca 



m R (A + wxa) Am R w a 2 Am R 
1 + log ^^ + ——^ log 



32vT 2 / 2 



A 2 



An 



327T 2 / 

m R w^ Nl a Am R w a 2 Am 

1 + ^ n log ~7~9 + 777T~o7o lo § 



'R 



A 2 

^AWI 



327T 2 / 2 

m R (A + w^ ecay 

2~?2 



16tt 2 / 



log 



327T 2 / 2 

Am R 



A 2 



A 2 

decay 



(48) 
(49) 
(50) 
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where 

, 10y 2 10 (m c V 2 \ 

w x = 6W + — -, w o = -p ( v ^— + (51) 

3 2 o 4 

A = ^l+M^l + 0(a). (53) 

Note that here m c /a = 0(1) and we recover the distinction among scale parameters (A, A , 

AawI Or Adecay)- 

These results reveal the following features of the WChPT. In general the chiral log 
terms(mR logm^) receive 0(a) scaling violation. In addition to this, the a 2 contribution 
generates logm^ term in m 2 , which is more singular as a function of m# than the usual 
chiral log term, rrtR logm^. Furthermore, both mjj logm^ and logm# terms are generated in 
fTT'Awi by the scaling violations, 0(a) for the former and 0(a 2 ) for the latter. The coefficient 
of logwiR term in uiawi is same as the one in m^. 

In the next section we employ the above formulae to fit the full QCD data obtained by 
the CP-PACS collaboration y. 



III. ANALYSIS OF CP-PACS DATA 

In this section, we apply the WChPT formulae to m 2 and mAwi i n the Nf = 2 full QCD 
with the clover quark action p]. 

A. Data sets and WChPT formulae 

The CP-PACS collaboration has performed the large scale full QCD simulations with 
the RG improved gauge action and Nf = 2 (tadpole improved) clover quark action, at 4 
different lattice spacings a and 4 different quark masses at each a, as summarized in table 
In ref. the data for m 2 and mAwi have been published. Unfortunately the data for f n at 
each quark mass are not available. 

We define the quark mass mjj in the WChPT theory in terms of the hopping parameter 
K in lattice QCD as 

Tfl Tfl I I 

m R = Z m (l + b m a—) — , ma = — -——, (54) 
m Uo 2K 2K C 



10 



where K c is the critical hopping parameter, and uq is the tadpole improvement factor, given 

( 8412 \ ^ 
1 — . This tur is identical to the renormalized VWI quark mass in 
P J 

ref. |8j. By definition, m 2 = at m R = in lattice QCD. We identify this m R in lattice 
QCD with wir in the WChPT, since m^, and therefore m 2 , must vanish at tjir = in the 
WChPT. We also use the renormalized wawi defined as 

m Awl = ^m%\. (55) 

2 



We employ the following fitting forms for m 2 and ttiawi 



m 2 = Am 



- 



R 



mRA + rriRawi ( AmR\ qwq [Am 



m A wi = A m R 



1 + 32n 2 f 2 bg V A 2 ) + 32tt 2 / 2 bg V A o 

rnjiawf^ 1 ( Am R \ a 2 w . / ArriR 
1 + ^7^ l0 HXi^J + 32^ l0g l^ 



IR 



(56) 
(57) 



B. Results 



We first fit the data at each a separately. Since there are only 4 data per observables 
at each a, it is impossible to fit an individual observable, m 2 or mAwi, as a function of rriR 
using eq. (j56j) or eq. ()57j) . each of which contains 4 or more parameters. Therefore, we try to 
fit m 2 and ttt-awt simultaneously. Since / can not be determined without data of /„., we fix 
/ = 93 MeV 15]. Even in the simultaneous fit, the number of independent fitting parameters 
is still too large. Since theoretically A Q = 1 in the continuum limit and the fit with A = 1 
becomes more stable, we fix Aq = 1 in our fit. In order to reduce a number of parameters 
further, we set Aawi = A = A, so we finally have 6 independent parameters, K c , A, A, wx, 
Wi WI and wo, for 8 data points. 

Fig. [T] shows data and fits for m 2 /mAwi as a function of ?tiawi a t eac h a. For comparison, 
the results by the fit with the standard chiral perturbation theory (wx = Wq = 0) are also 
given. It is manifest that the WChPT fits perform much better than the ChPT fits. The 
parameters extracted from the fits are given in table |H] Note however that xVdof shown 
in the table has not been reliably estimated due to the correlation between m 2 and mAwi, 
which is not given in ref. [sj]. 

In Fig. El A, A, Wxa, Wi WI a and w a 2 are plotted as a function of a, together with K c as 
a function of the bare gauge coupling constant g 2 . While A, A and W\a are too scattered to 
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be fitted, K c , w$a and wf a may be fitted as 

K 1 1 + d (K c )g 2 + d^K^ + 4(^ c )^ 6 , s 

8 ' 1 + (do(K c ) - 0.02945)# 2 1 j 

where 0.02945 is the 1-loop coefficient and 

wf WI a = do(wx)a, w^a 2 = do(wo)a 2 . (59) 

Fit curves are also shown in Fig. El and the extracted parameters are given in the column 
(a) of table EH 

To determine a dependences of A, A and w±a, we have fitted m 2 /mAwi as a function of 



both m R and a, using the following formula derived from eqs. (j56|57j) with Aawi = A: 



(60) 



(A + Awia)m R ( Am R 



32tt 2 / 2 V A 2 



where 



A = d (A) (I + d 1 (A)a + d 2 (A)a 2 j , A = 1 + d Q (A Q )a (61) 
A = do(A) (l + rfi(A)a 2 ) , A Wl =w x - wf l = d {A Wl )a. (62) 



No \ogm R term is presented in eq.(|5D|l. Note however that logm-Awi term appears again 
if we replace m R in the right-hand side of eq. lpjUj) with wawi, due to the presence of the 
\ogm R term in eq.([57|). With K c fixed to the values in tabled the fit works well, as shown 
in Fig. |2l and the fitted parameters are given in the column (b) of table 11111 

We roughly estimate the size of each parameter, Bi, V 2 , Wj.,2,3 from the continuum 
extrapolations of A, w\, wo and m c . Since we can not separate the 1/a contribution in 1/K C , 
however, m c can not be extracted. Therefore, we simply set m c = 0, giving that W% = 2aW 2 ; 
the leading contribution of W\ vanishes. To reduce the number of the parameters further, 
we set Wo = 0. Then extracting B\, W 2 and V 2 as 

B, = f 2 d (A) = (A Bl ) 3 (63) 
W 2 = s { A W2 f (64) 

_ f 2 d ( Wl ) _ f 2 (d (wt wl ) + do(A Wl )) 4 
2 10 10 1 ' 

we obtain A Bl = 0.41 GeV, A Wa = 0.24 GeV and A V2 = 0.21 GeV. These A x takes a 

reasonable value, Ax = 0.2 ~ 0.4 GeV. If aAx > m/Ax, O(a) terms become more important 

than m R terms. With A x = 0.2 ~ 0.4 GeV, this condition at a~ l = 1 GeV or a _1 = 2 GeV 

corresponds to m R < 40 ~ 160 MeV or m R < 20 ~ 80 MeV, respectively. 
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C. Validity of the (W)ChPT 



We now estimate the relative size size of the next-to-leading contribution to the leading 
contribution in the WChPT for m%: 

n/ wp ,nrr\ m n( A + oWj) + a 2 w ( Am R \ 

fl(WChPT) = log ) (66) 

for the WChPT at finite a, where parameters A, A, w± and Wq depend on a. We plot 
i?(WChPT) in Fig. Hat a(GeV- 1 ) =0, 0.44(0 = 2.2), 0.55(0 = 2.1), 0.79(0 = 1.95) and 
1.1(0 = 1.8). While the 1-loop contribution takes reasonable values, 10% ~ 30 % , at 0.1 
GeV < m R < 0.2 GeV for all a, the contribution from logm^ in the WChPT diverges 
as m R — > 0. This might invalidate the WChPT in the chiral limit. We will consider this 
problem in the next section. 



IV. RESUMMATION OF log m R TERMS 

As evident from the analysis in the previous subsection, \ogm R contribution becomes 
larger and larger toward the chiral limit, so that we can not neglect "higher order" term 
such as (logm^)™ (n = 2,3, •••). We must perform a resummation of logm^j term at 
all orders. Since it is possible in principle but difficult in practice to calculate (log m R ) n 
contribution at n-loop order jWe derive resummed formulae from a different point of view. 

As discussed in refs. IjJ, the massless pion corresponds to the inverse of the divergent 
correlation length at the second order phase transition point. Since the effective theory 
which describes this phase transition is some 4 dimensional scalar (pion) theory with rather 
complicated interactions jl^] . the phase transition has the mean-field critical exponent with 
possible log-corrections. In particular the pion mass, the inverse of the correlation length, 
should behaves near the critical point as 



rat 



Cm R [log +•••, (67) 

where • • • represent less singular contributions. If we expand 



+ + (68) 
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where 



A' = log (§) (69) 

the formula at the next-to- leading order in WChPT, eq.([56jh is recovered, with the identifi- 
cation that 

v' _ a 2 w , _ 

X ~ 3W " ( } 

To determine v' and X separately, the explicit calculation in the WChPT at 2-loop or more 
orders is necessary. This will be considered in future investigations. 

We have finally obtained the following resummed formulae for m\ and tuawi'- 



ml = Am R \\og (^)} 32 " 1 



aP'w, 



a wq 

m Awl = A m R jlog ("J^)} 32 " 1 



m R A + m R awi ( Am R 
1 + 327r 2 / 2 lQg l^ 



mnawf^ 1 ( Am R 



32tt 2 / 2 V A a 



AWI, 



(72) 
(73) 



where A, A and ojq may be different from those in eqs. f|56|57j) . It is better to use these 
formulae instead of the previous ones, eqs. ()5fi|57j) . in future investigations. Eq. (|fiUJ) remains 
the same. 

As a trial, we use these formulae with Aq = 1, Aawi — A and A =l GeV, in order to fit 
ml and mAwi simultaneously, at each a. The quality of the fit is as good as the previous 
one, and the fitting parameters are compiled in the end of table |H] In addition, the next- 
to- leading contribution, the second term in eq. (|73jl . vanishes toward m R = as shown in 
Fig. HI where i?(WChPT) in the previous subsection, which is now modified as 

n,,„n, DT v m R A + m R a Wl ( Am R \ 

i?(WChPT, resum) = log ) , (74) 

are plotted at = 1.8, 1.95, 2.1 and 2.2 . 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper we have derived the effective chiral Lagrangian which includes the a 2 effect 
of the Wilson-type quark action in the case of the Nf = 2 degenerate quarks. Using this 
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effective Lagrangian the quark mass(mfi) dependences of m%, ^awi and f n have been cal- 
culated at the 1-loop level. We then have simultaneously fitted and ttiawl, obtained by 
the CP-PACS collaboration for Nf = 2 full QCD simulations, using the WChPT formula, 
and have found that the data are consistently described. We have attempted the continuum 
extrapolation of the WChPT formula. 

Comparing to the standard ChPT, several distinct features such as the additive mass 
renormalization, 0(a) corrections to the chiral \og(mn log tjir) term, a more singular 
term(logm#) generated by 0(a 2 ) contributions and the presence of both m^logmR and 
logm# terms in itiawi, leads to the success for the WChPT formula to describe the CP- 
PACS data. Although an ambiguity for the definition of K c caused by the additive mass 
renormalization can be avoided by the use of mawi, the last feature, the existence of both 
rriR logm^ and logm^ terms in mAwi, makes the WChPT formula different from the ChPT's. 
The large 0(a) correction to m^logm^ term plays an essential role to describe the actual 
data, though more or less others have some contributions. We have also derived the formula 
after resumming log ttir terms, using the fact that the mean-field critical exponent receives 
the log- correct ion. 

Because of the limitation of available data, our WChPT analysis is far from complete. 
Therefore it is important to refine the analysis by taking the correlation between m 2 and 
mAwi m t° account and including f n data in the simultaneous fit, in order to establish the 
validity of the WChPT. Reanalyses of other full QCD data have to be done of course. It is 
also urgent to derive the WChPT formula for other cases jl^ such as the quench/partially 
quench cases, the Nf = 3 non-degenerate case, vector mesons and baryons, heavy-light 
mesons. 

Once the validity of the WChPT to describe lattice QCD data is established, instead of 
thinking that the quark masses in the current full QCD simulations are too heavy for the 
ChPT to apply, we may say that some (but not all) of lattice data are well described by the 
(Wilson) chiral perturbation theory, by which errors associated with the chiral extrapolation 
may be well controlled [l^j ]. 
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TABLE I: Parameters of N f = 2 full QCD simulations by the CP-PACS collaboration 8]. The 
scale a is fixed by m p = 768.4 MeV. 



P L 3 x T csw a [f m ] a 1 [GeV] La [fm] m^/nip 

1.80 12 3 x 24 1.60 0.2150(22) 0.9178(94) 2.580(26) 0.55~ 0.81 

1.95 16 3 x 32 1.53 0.1555(17) 1.269(14) 2.489(27) 0.58~ 0.80 

2.10 24 3 x 48 1.47 0.1076(13) 1.834(22 ) 2.583(31) 0.58~ 0.81 

2.20 24 3 x 48 1.44 0.0865(33) 2.281(87) 2.076(79) 0.63~ 0.80 



TABLE II: Parameters of the WChPT fit at each (3. 



p 


K c A [GeV] A [GeV] Wl a[Ge\] w a 2 [GeV 2 ] w^ l a [GeV] 


rVdof 


1.80 
1.95 
2.10 
2.20 


0.147761(15) 5.114(28) 0.079(19) -5.525(64) 0.206(22) -0.560(74) 
0.142160(19) 5.377(33) 0.193(51) -5.162(74) 0.241(42) -0.457(118) 
0.139110(12) 5.807(14) 0.694(20) -5.24(18) 0.417(50) -1.15(27) 
0.137691(23) 5.669(71) 0.128(88) -5.15(20) 0.039(16) -0.22(39) 


0.3 
0.3 
0.2 
0.7 




resummed WChPT 




1.8 
1.95 
2.1 
2.2 


0.147562(15) 5.111(29) 0.067(12) -4.862(46) 0.787(21) 0.124(15) 
0.142009(7) 5.366(23) 0.132(15) -4.538(52) 0.624(18) 0.310(32) 
0.138959(13) 5.535(47) 0.131(71) -4.79(14) 0.280(37) 0.181(49) 
0.137657(36) 5.789(106) 0.391(82) -4.63(12) 0.201(95) 0.195(96) 


1.5 
0.3 
1.2 
0.8 
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TABLE III: Continuum extrapolation of the WChPT fit parameters, (a) and toawi are fitted 
as a function of itir at each a. Then parameters are fitted as a function of a. (b) m^/rrtAwi are 
fitted as a function of itir and a 



(a) 


(b) x7cfof=1.3 




d (X) d ± (X) d 2 {X) x 2 /dof 


X 


d (X) d!(X) d 2 (X) 


K c 
w 

w AWl 


-0.2127(10) -0.008300(55) 0.000787(31) 3.6 
0.202(17) 20 
-0.549(61) 3.4 


A 
A 

^0 


8.087(97) -1.002(29) 0.2672(29) 
1.196(35) -0.8404(58) 
-1.62(25) 
-0.590(47) 
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FIG. 1: The WChPT fits for and toawi a t eacn P- Results are shown for to^/toawi as a 
function of tbawi- For comparison the standard ChPT fits (wi = wq = 0) are also included. 
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FIG. 3: The WChPT fits for m^/mAwi as a function of uir and a. Results are shown for m^/mAwi 
as a function of uir. 
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FIG. 4: Left: The relative size of the next-to-leading contribution to the leading one in the 
WChPT as a function of the quark mass tur at /3=1.8,1.95,2.1 and 2.2, together with the one in 
the continuum limit (ChPT). Right: Same quantities in the resummed WChPT. 
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